Relativistic electrons in a rotating spherical magnetic dipole: localized 

three-dimensional states 
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Paralleling a previous paper, we examine single- and many- 
body states of relativistic electrons in an intense, rotating 
magnetic dipole field. Single-body orbitals are derived semi- 
classically and then applied to the many-body case via the 
Thomas-Fermi approximation. The many-body case is rem- 
iniscent of the quantum Hall state. Electrons in a realistic 
neutron star crust are considered with both fixed density pro- 
files and constant Fermi energy. In the first case, applicable 
to young neutron star crusts, the varying magnetic field and 
relativistic Coriolis correction lead to a varying Fermi energy 
and macroscopic currents. In the second, relevant to older 
crusts, the electron density is redistributed by the magnetic 
field. 
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1. INTRODUCTION 



In a previous paper Q (hereafter paper I), we exam- 
ined the relativistic semiclassical relativistic orbitals of a 
particle of mass m and charge q in an intense magnetic 
field, idealized as a dipole, in a rotating reference frame. 
The particle was confined to the spherical surface. In this 
paper, we present the treatment of three-dimensional or- 
bitals, the local cyclotron or Landau states. These re- 
sults are applied to many-body electron states defined in 
a semiclassical (Thomas-Fermi) approximation and cal- 
culated in a simplified neutron star crust model with elec- 
trons and nuclei. We consider quantum dynamics only 
for the electrons, not the hadrons ||^, but effects of rela- 
tivity and rotation are included. 

We ignore gravity (derivatives of the metric) here, as 
this is negligible for charged particles compared to the 
magnetic field. Where numerical values are needed, the 
tilt angle between the dipole and rotation axes is as- 
sumed to be maximal, sin^o Ij and the rotational ve- 
locity to be w = VlR/c — 0.01, a realistic value for 
high- field neutron stars with radius R . 

The treatment is based on expanding the particle mo- 
tion in inverse powers of the field strength. Although 
electrons are stripped from the neutron star surface by 



the rotation-induced electric field, the bulk of electrons 
remain in the crust to preserve local charge neutrality, 
with the surface sheathed by a thin space charge. The 
space charge is stabilized by the Coulomb force (with 
the positive crystal) opposing the induced electric field. 
Only a small fraction of electrons are accelerated into 
the stellar wind Our final applications here are to 
quantum single- and many-body states where radiation 
emission is neglected. This is exact for charged parti- 
cles in their ground states or in excited states unable 
to decay by Pauli exclusion blocking in the presence of 
other fermions. We also seek a general classification of 
possible orbitals based on the relevant kinematic param- 
eters. The semiclassical quantization is based on the 
Wilson- Sommerfeld or Bohr-Sommerfeld rule, a result of 
the WKB approximation p]. 



FIG. 1. Geometry of the magnetic dipole M and sphere 
rotating at angular velocity fJ, with relative tilt Oq. 



2. GENERAL RELATIONS 

2.1 Coordinates, metric, and field 

As the magnetic field dominates rotational effects (un- 
less r is quite large), we place the rotation axis at an 
angle do with respect to the magnetic dipole in the cf) — 
0, TT plane (Fig. 1). The metric in a spherical polar coor- 
dinate system (r, 6, (p) rotating with the sphere is given 
by the line element 
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ds = g^ijdx^dx" 

= c2(l - u}^)dt^ - dr^ - r^{de^ + sin^ 
2crujrj, sin dt d(j) — 2crujg dt dO 



(1) 



The vector u is defined from the rotational angular ve- 
locity vector r2 by a; = X r/c. We also use uj = ilR/c, 
where r = i? is a reference sphere. The components of 
CO are 



W0 = {ujr/R) [cos 9o sin 6 — sin cos 9 cos ( 
We = — (Ljr/i?) sin^o sini/) , 



(2) 
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These are the appropriate generalizations of paper I and 
the treatment of Landau & Lifshitz Q to the case sin Oq ^ 
and r ^ constant. The lightsphere is the surface CJ^ = 
1, or 



up-{rlRf\cos^ 6*0 sin^ 0+ 

sin^ 00 sin^ ^ sin^ 0o cos^ Q cos^ ( 

(1/2) sin 26*0 sin 26* COS./)] = 1 



(3) 



This is the surface upon which 500 = (Fig. 2). 

We choose our axes so that the magnetic dipole is along 
the = direction. The dipole magnetic field has polar 
strength at r = i? : 



v4e = 
= 

2r 



(4) 
(5) 

(6) 



where is covariant in the rotating spherical coordi- 
nates. For convenience, we rescale Bq into dimensionless 
form as 



/3o - k|Soi?/(2mc2) . 
The magnetic moment of this dipole field is 
|M| = . 



(7) 



(8) 



In a high-field neutron star, /3o ~ 10^^. 

The Lagrangian of the charged particle is expressed in 
terms of the proper time r, with = dx'^ /dr. There are 
four equations of motion including one for each momen- 
tum component and one energy-momentum constraint. 
Since the Lagrangian does not explicitly depend on t, we 
have 



dL 
'dt 



= 



which implies 



Because 



dPo^ 

dr 



(9) 



(10) 



dr = 

we have also 



dt\/9t^y {dx^' /dt) {dx" /dt) 



dPo 
dt 







(11) 



(12) 



where Pq = E, the energy. The equations of motion for 
r, 0, and 9 are non-trivial: 



dL dPr dL dP,_ 



dL dPg 
'd9^~d^ 







(13) 



dr dr d(j) dr 
The energy-momentum-mass constraint is 

g°0p2 + 2g°'Po(P. - ^A,) + g^^iP, - ^A,)(P, - ^A,) 

= (mc)2 , (14) 

including radial and angular terms, and the mixed rota- 
tional g^^ , g^'^ terms. The contravariant metric compo- 
nents are 

.9°" - l/c^ , r = -1/r^ , 

= -(1 - ^l)/r^ , 9^^ = -(1 - ^l)/{r^ sin^ 9) , (15) 



-W0/(C 
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'Uj0/{cr) . 



Along the actual worldpath in spacetime the energy- 
momentum constraint is equivalent to g^^ xt^x'' = c^. This 
condition is valid after varying the action and simplifies 
the equations of motion. 



2.2 Asymptotic orbits 

As a check of this dynamical system, we consider 
briefly the r ^ 00 orbits. The constraint (14) alone is 
sufficient to give essential information about the trajec- 
tory. Note first that uP' < 1 for physical trajectories. 
Since ui^ — ^"a, + each component must also have 



are never 



magnitude less than one. Then g^^ and g'f 
exactly zero. 

In paper I the natural energy scale for unconfined 
charged particles at the surface r = R was set by e ^ /3o, 
or E ^ \q\BoR/2, independent of m (see also Longair ^). 
In the limit r 00, Eq. (14) simplifies to 



p2 



(E/c)^ ~ {mc)' 



(16) 



where E is naturally <; 0{\q\BoR). In this regime, Bq 
is much smaller than at the surface of the star, as r is 
much larger. The magnetic potential term vanishes as 
l/r^, while the angular momentum terms vanish as 1/r. 
Particles can escape as r ~> cxd if > and if they 
remain within the lightsphere surface ui^ = 1 (Fig. 2). 
Both graphically and analytically, it is seen that {(j>, 9) 
(O,0o) or (tt, TT — 9o), as r 00. That is, an escaping 
particle is forced into the common plane of M and f2 
and leaves along the rotation axis. 
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FIG. 2. Lightsphere r{9,(f>) defined by Eq.(3), surface of 
goo ~ 0, for cu — 0.01. 6 and (j> in radians. 

For sufficiently large r (strictly speaking, for suffi- 
ciently large | the magnetic dipole approximation 
breaks down, and the field (6) can no longer be used. 
The field lines are twisted about by the rotation near 
the lightsphere, and the particle motion is more com- 
plex ||9(]. Gravitational attraction is negligible compared 
to the magnetic force for such energetic particles, even 
far from the sphere. 



3. LOCAL SINGLE-BODY STATES 

The local Landau states can be obtained by expand- 
ing the energy-momentum constraint in a Taylor series 
around r = r, i.e., let r go to f -|- Ar. The local point 
is defined by its spherical coordinates (f, 9, (j)) . We retain 
terms up to second order in Ar and momentum P, but 
leave the local metric g^i,y{v) constant. 



3.1 Local coordinates and field 

We define a local Cartesian coordinate system so that 
the z co-ordinate is along the local magnetic field vector, 
the X coordinate is orthogonal to z and pointing outward, 
and the y coordinate is azimuthal. 

The dipole magnetic field is 



B = (So/2)(-) [2 cos 6* f + sin 61 d] 



(17) 



In terms of the local displacement Ar, the new coordi- 
nates are given by 

Ar sin 9 - 2rM cos 9 



Vl + 3cos2 6l 
2 Ar cos 9 -f r A^ sin 9 



\/l + 3cos2 6l 
y — — r sin 9/S.<f> , 

while the local Cartesian basis is 

sin 9 Y — 2 cos 9 6 



(18) 



z = 



y/l + i C0S2 9 

2 cos 6* f + sin 9 d 



VI + 3cos2 6i 
y = -"^ . (19) 
The bar on r is now omitted unless needed. 

3.2 Local Landau states 

The constraint can be rewritten in the local coordi- 
nates as 



e^-nl-nl- (n; - f3x/rf - 1- 



2uJ4,€{n' - (3x/r) - 2ujge 



2ujAe\ 



j3 cos ( 



sin - 2 cos 6* W^, 
Vl + 3cos2 6l 

[4 cos 9 x"^ - 



\r2sin6l(l-H3cos2 6l)3/2 
3 cos 9 sin^ 9 z'^ + 2 sin 61(3 cos^ 6* - 1) a;z] 



(n; - px/rfiol + 

(sin 6* - 2 cos 61 Ha; ) 2 2 



1-1-3 cos2 9 



(20) 



Here we use dimensionless local magnetic field strength, 
energy, and momenta given by 



2mc2 

e = E /mc? 

^x,z = Px,z/mc 



(21) 



n' is the new canonical momentum given by the trans- 



formation 



n' 



qB^R^ sin 9 



(22) 



The second term is a shift which is constant in the local 
coordinate system. The first line of the constraint (20) 
gives the classical version of the Landau system in terms 
of the local magnetic field B(r), which also defines the 
local z axis. The other terms are linear and quadratic 
corrections (Coriolis and centrifugal terms) due to the 
local u components uj^ and loq. 

Neglecting rotational effects yields the classical cy- 
clotron motion [Q. In the semiclassical form, these are 
the Landau orbitals, with energy eigenvalues 



i + nf 



— ) (^) [2ni + 1 - sgn(g)a] , (23) 



where the principal Landau quantum number is n^, = 0, 
1, 2, ... and the spin a has been included, cr = ±1 P,p^. 
For sgn{q) < 0, the state {nL,cr = 1) is degenerate with 
{riL + 1,(7 = —1). The Landau orbitals are character- 
ized by length and momentum scales L = \J'hcl (2|gB|) 
and Pj^ — y^h]qG\/{2c) in the {x, y) plane transverse to 
the field. The longitudinal momentum — mcll^ is an 
eigenvalue continuous over the range —oo<Pz< +oo, as 
the motion in that direction is force-free. Note that the 
magnetic factor proportional to f3 is suppressed by the ra- 
tio of the Compton wavelength Ti/{mc) to the spherical 
distance r. The effect of the magnetic field in the en- 
ergy is thus negligible compared to mc^ unless (3 is very 
large, as seen in paper I. Single-body Landau states of 
fermions are unstable against radiation, unless the lower 
energy states are already filled, as they are in Sect. IV 
below. 

We now consider the effect of rotation as a first-order 
perturbation in (D, approaching the problem semiclas- 
sically. (The terms quadratic in uj are much smaller.) 
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The unperturbed orbitals are Landau states in the {x, y) 
transverse plane. We keep the terms hnear in ui (Cori- 
oHs effect) and average them over one Landau orbit, as 
in paper I, Sect. V. But unUkc that case, this averaging 
is over a local microscopic orbit, not a macroscopic or- 
bit over the whole sphere. The terms linear in u) thus 
do not average to zero in general. The modified classical 
energy-momentum constraint is then 



(n' - px/rf 



l- 



2ujge 



sin 9 n. 



VI + 3cos2 6i 

/3cos( 



r2 sin 61(1 -h 3 cos2 61)3/2. 
4cos6l - 3cos6'sin2 6' z'^] = 



(24) 



noting that {x)o = {y)o = (nx)o = (ny)o = 0, averaged 
over an unperturbed cyclotron orbit. These orbits have 
semiclassical radius a{nL) = 2L^nL + 1/2. The semi- 
classical result for the Landau energies corrected through 
0(u}) is then 



2cjeeosin6' n^/vl + 3cos2 



Vr^s] 



/5 cos ( 



^sm«(l + 3cos2 61)3/2 
[4(2nL -f 1) cos 61 - 3 cos 6* sin^ 61 



cJ 



(25) 



where we have used (112)0 = H^, (a;2)o = a?{nL)/2, and 
a distance cutoff zc for motion in the z direction. (This 
cutoff is discussed further in Sect. IV below.) The last 
term can be neglected if the Landau orbits are much 
smaller than the spherical distance r : L, zc « as 
it is another factor of L/r smaller than the zeroth-order 
term. The energies depend on 112 zeroth order but 
receive a correction proportional to II^ at first order in 
CJ, breaking the symmetry — > — H^. The Landau pole 
states of paper I can be recovered if sin 6' ^ 0, r = i?, 
identifying with the old ng and setting = 0. 

The error arising from neglecting the cubic terms in 
the Taylor expansion can be estimated and varies as 
~ |Ar|/r times the quadratic terms. That is, the cu- 
bic and higher terms in the expansion are suppressed by 
additional powers of the Landau length L over the sphere 
size r. 



4. LOCAL MANY-BODY STATES 

4.1 Density of states 

The full density of states is a product of four factors: 
the density of Landau states, the degeneracy factor 
in the (a;, y) transverse plane for each Landau state, the 
density of longitudinal states for z motion, and the 
spin factor (one for the ground state — Q, two other- 
wise) . 



The degeneracy of a given Landau state in the 
transverse plane is 



|gB| 
2ttTlc 



(26) 



per unit planar area, a result valid in both non- 
relativistic ||l^ and relativistic regimes, while the lon- 
gitudinal motion contributes a factor 



Vn 



nic dHz 
2Trh 



(27) 



per unit longitudinal length. Thus the number of states, 
including the spin factor, is 



dS dz 



{2-5r,^fl)-V^-V. 



(28) 



per unit transverse surface area dS and unit longitudinal 
length dz. 

In the semiclassical limit, where n^, is quasi- 
continuous, the density of Landau states per unit energy 
is given in dimensionless form by 



dn 



[1 



i finer 

dT ~ yJpn 



) [1 + (2cje/v/l + 3cos2 9)]- 
2hf3/{mcr)\^nL -f 1 — sgn{q)a] 



2bJetQsin9 [Hz/Vl + 3cos2 



(29) 



4.2 Thomas- Fermi approximation 

The Thomas-Fermi method approximates quantum 
many-body fermion states in a varying potential with 
local states defined by a locally constant field Q| . In our 
case, the local electron states are filled up to some highest 
and partially-filled Landau level nj^ by the charge carri- 
ers, assumed here to be electrons. That part C,"^ of the 
squared energy arising from the 11^ terms alone, 



C'(n.) = n2 



2uJa 



1 + 3 cos2 9 



must be cut off at some maximum, as must the corre- 
sponding z motion. In a real system, the cutoffs are 
provided naturally by the presence of lattice ions: |zc| ^ 
h/ {Zcsctrac) (Bohr length) and jllzl ^ Z^sa, where Z^s 
is an effective (screened) positive ionic charge |^ . In the 
realistic case the zq term is unnecessary; assume that 
< Cf = (ZcSOi)^. Then the longitudinal momentum 
Hz is cut off asymmetrically at Hz = 11+ > and < 
0, with n+ ^ -n- if 7^ (Fig. 3). 
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FIG. 3. Momentum-dependent part C,"^ of the squared en- 
ergy as a function of longitudinal momentum Hz , including 
0{ijj) correction, showing the solutions 11 J of = Cf- 



by 



The electron number density rig at any point is given 



kB| 

2'Khc 



[l + 2{nl-l) + 2iy]- An,imc/2Trh) , (30) 



where is the partial filling factor of the highest Landau 
level n*]^, < I' < 1, and is assumed > 1. (If = 
0, the entire term in brackets is replaced simply by v.) 
Because of the z degree of freedom, each Landau level is 
actually a band, with lowest energy level at some Hz ^ 
0, whose sign depends on that of ojg. In addition, there 
is the planar degeneracy, modified by the partial filling 
factor v in the highest Landau level. If rig is specified 
along with AII^, then nj^, i/, and ep can be determined. 
Procedural details are found in the Appendix. 

At zero temperature (the case examined in this paper), 
the electrons cannot radiate into already-filled lower- 
energy single-body states. Radiation occurs only if 
fermions are externally excited above the Fermi energy. 

If He is held constant over a sphere of radius r, then 
n*j^, v, and ep change over that surface as |B| changes. 
The variation of and over the surface is given by 



in turn depends on the nuclear crystal and the mag- 
netic field. (Relativistic electrons travel at essentially 
the speed of light, leading to saturated current densities 
of ©(riec), apart from magnetic effects on the density of 
states As the magnetic field and rotation affect the 
electron Fermi energy, the beta equilibria of neutrons IT^ ] 
and muons are also affected above their respective thresh- 
olds. This last effect cannot be included without revising 
the hadronic equation of state . The Fermi energy will 
be unequal over the sphere at first, a limit expected to 
apply to young neutron star crusts. With time, ep will 
equilibrate to the same value everywhere, and the B and 
He, rip profiles will change in parallel. 



Int 



,2Tih 



V ^r. I V ATT / V 



mc 
1 



- 3 cos^ ( 
27r7ix3/ rig 
mc ' ' Alia 
1 

4Vl + 3 cos2 ( 



Vatt I^Bo^' 

1 * 
rii 

2 ^ 



(31) 



where Int[/] in (31) denotes the integer value of the func- 
tion /. Typical variations of and v over the sphere are 
shown in Figs. 4 and 5. The Fermi energy at any point 
on the surface is 



i+ep + {^ 

Vl + 3cos2 6' 



f ) {2nl + l + a)x 



(32) 



where sgn(q) = — 1, and ap = —1 for the level; the 

critical field strength is Be = [ml / n) ^J^FJh^ = 1.3 x lO^^ 
G, defined in paper I. Matching variations of ep are 
shown in Fig. 6. As the Fermi energy ep varies over 
a sphere, electron currents flow as the electrons seek 
the lowest energy. There are also radial currents (see 
Sect. IV.C below). 

Quantitative details of the subsequent evolution p^ ] 
are beyond the scope of this paper, but certain qualita- 
tive features are clear. The magnitude and evolution of 
the currents depend on the electrical conductivity, which 



FIG. 4. The last, partially-filled Landau level on a 
spherical surface, at constant density, defined by Eq. (31); Bo 
= WBc. (a) p = 10 g/cm^; (b) p = 10^ g/cm^; (c) p = lO" 
g/cm^. 9 and (f> in radians. 
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FIG. 5. The partial filling factor v of the last Landau 
level nj, on a spherical surface at constant density, defined 

by Eq. (31); Bo = WBc (a) p = 10 g/cm^ (b) p = lO'^ 
g/cm^; (c) p = 10^^ g/cm^. 9 and 4> in radians. 



FIG. 7. The electron number density rie on spherical sur- 
face at constant Fermi energy, defined by Eq. (30); Bo = Be. 

Units arc inverse Compton volume, (a) ep = 1.06; (b) €f = 
5.0; (c) = 1800. 9 and 4> in radians. 



4.3 Radial structure 



For a simple neutron star, we can assume a given profile 
of positive ions, and thus electrons, with radius r. With 
the field |B(r)| in addition, the radial dependence of ep 
can be found. 

In Fig. 8 the radial profile is shown for ep from r/R = 
0.7 to 1.0, with an expanded subfigurc oi r/R — 0.9998- 
1.0. Even at this detail, the thin surface non-relativistic 
layer cannot be seen, but the quantum regime of dis- 
cretizcd ep steps is clearly visible. Fig. 9 shows ep as a 
function of r/R and 9. The rotational correction drops 
out in certain cases. If sin</) = (r in the M-fi plane) 
or sin^o = (no tilt), toe vanishes. Also, if sin^ = 
(f along the dipole axis), Pq ^ rsinO ■ 0, and the 

Coriolis effect disappears. Otherwise the ep radial profile 
is affected dramatically by the rotational correction for 
relativistic ep. 



FIG. 6. The dimensionless Fermi energy ep on spherical a 

surface at constant density, defined by Eq. (32); Bo = lOBc- 
(a) p = 10 g/cm^ (b) p = 10^ g/cm^; (c) p = 10^=^ g/cm^ 9 
and (j) in radians. 

The opposite limit is the case where the Fermi energy is 
constant in space and no currents flow, a situation prob- 
ably holding for crusts at later times. Since the field |B| 
varies, the electron density rig, and thus the positive ion 
density n^, must also vary. This is the case where the 
magnetic field is strong enough to dominate the mechan- 
ical structure of dense matter. We show three examples 
of how rie is redistributed over a sphere of constant r for 
fixed ep (Fig. 7). The proton density is rip = (Z/Zgfj)ne, 
and p can be inferred from Eq. (A2). 



FIG. 8. Dimensionless Fermi energy ep as function of r/R 
for simplified neutron star crust model of text; sin^ = sincf) = 
1, and Bo = WBc. (a) r/R = 0.7-1.00; (b) r/R = 0.9998-1.00. 
Terraced steps indicate quantum regime at the surface. 
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FIG. 9. Dimensionless Fermi energy ep as a function of 
r/R and 6 for simplified neutron star crust model of text, (a) 
sin0 = 1, r/R = 0.7-1.00, Bo = lOBc, smO > 0; (b) detail 
showing sin6' > 0.2. 9 and <j) in radians. 

For this crust profile, we have used a neutron star 
model of Glendenning with only electrons, nuclei, 
and neutrons, and no muons or other hadrons. The ra- 
dius and mass of the neutron star are i? = 11 km and M 
= 1.55M0, respectively. The profile starts at r/R = 0.7, 
where p — 2 x 10^^ g/cm^ and decreases outward. As 
tF varies in space, currents flow in response, a scenario 
probably relevant to early neutron star crusts. (Recall 
that we ignore the effect of gravity.) Again, the neutron 
and muon beta equilibria are changed from their zero- 
field analogues. 

Now consider the opposite limit again, probably hold- 
ing for older crusts, a constant ep, with ne{r) determined 
by a given field profile B(r). Then v must be set inde- 
pendently; here we take v = 0.1 and 0.9 as illustrative. 
For critical field strength, the density rises with r, while 
for even higher fields, the density is quenched to essen- 
tially a constant (Fig. 10). For the same surface value of 
e_F, the Tie density of Glendenning ||^ is too large to be 
shown on the scale of Fig. 10. Thus, the crust densities 
are low compared to the non-magnetized values when the 
magnetic field is the controlling factor. 



with chemical equilibrium and Coulomb neutrality. De- 
pending on the degree of lattice disorder and interelec- 
tron forces, various conducting, insulating, or quantum 
Hall-like many-body states can arise [|ll|,|lj], affecting the 
macroscopic currents inferred in Sect. IV. 
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APPENDIX: MANY-BODY THEORY 

To determine the filling of momentum states labelled 
by n^, the cutoff equation (^"^{Jlz) = Cf niust be solved 
for some given value of (p, taken here as (Zcfja)^. At 
zeroth order in (D, this equation is trivial. Once the rota- 
tional corrections are included, the equation is not only 
quadratic and linear in 11^, but has an implicit depen- 
dence on Hz through eo- 

In the 0{u!) correction, our procedure is to take eg = 
1 + where e\ is the purely two-dimensional Landau 
term: 



1 - sgn(g)cr] 



FIG. 10. The electron number density as a function of 
r/R for constant ep = 1.06; sin = sin (j> — 1, and Bo = lOBc. 
Solid: V — 0.1; dotted: v — 0.9. Density profile of simplified 
neutron star crust model in text is too large to show on this 
scale. Units are inverse Compton volume. 



5. CONCLUSION 

This concludes the treatment of semiclassical orbitals 
begun in paper I. In this paper, we have found the local 
one- and many-body states of relativistic charged parti- 
cles confined to a sphere with an intense, rotating mag- 
netic dipole field. 

There remain full quantization with the Dirac equation 
and the inclusion of the positively-charged lattice struc- 
ture to determine the local neutron star matter state. 
A full calculation requires a self-consistent treatment of 
gravity, nuclear matter, magnetic fields, and currents. 



and neglect the 11^ term, as the latter is typically much 
smaller than one. In that case, the cutoff equation ('^{Hz) 
= is a simple quadratic with the two roots and 



An, 



^^ 



n: 



2a/C| + + el) sin^ 9/{l + Scos^ ( 



(33) 



as the allowed spread of z momenta. Given a value of rie, 
the values oi ep, n*^, and v are found iteratively, starting 
with u) — Q, then with this solution used in the 0{u)) 
corrections. 

In the opposite case, of fixed ep, n*j^ is determined, 
while we set i' = 0.1 and 0.9 as illustrative (Fig. 10). 
The density rig is then found. As long as v ^ 0.5, 
can fall or rise with r. A complete treatment of the bulk 
crust requires inclusion of nuclear matter and gravity [D , 
as well as an interior B field profile, not necessarily a 
dipole. 

A semi-realistic spatial profile of electron density re- 
quires the proton density Up = n^, usually determined in 
terms of mass density p. The "effective" electron density 
Tie, the density available for conduction, is 
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(2TTh mc) ne^cond = — j ^ — , (34) 

A g cm"'^ 

where Z^s is the number of electrons per nucleus available 
for conduction (unbound electrons) . The atomic number 
A = N + Z, where the neutron number N per nucleus 
is abnormally large for nuclei in an electron Fermi sea. 
The number density is normalized to a Compton volume 
{2TTh/mc)^. From p = 10 to about 3 x lO'^ g/cm^, the 
inner electrons of the atoms remain bound, not partici- 
pating in conduction; in this case, Z^g < Z and can be 
read off from standard atomic structure jist . The nu- 
clei are always iron {Z — 26 and Z^s = 8, 16, 24) at 
these densities P,p^. For higher densities, the orbitals 
of different atoms merge, and Z^g = Z. (We neglect the 
formation of partially ionized atoms at ultrahigh densi- 
ties, assuming that all electrons are stripped from their 
nuclei.) The results do not depend sensitively on Zcs- 

For nuclear matter composition at densities below com- 
plete nuclear dissociation, but for Z > 26, we use the 
results from |l|,0 (see also ^) at p < 2 x 10" g/cm^. 
Our method can be applied at higher densities with 
free nucleons, but it needs to incorporate the presence 
of muons and then of heavier strange and non-strange 
hadrons, and then possibly of a quark-gluon plasma pi. 
In this paper, only the simplest case of electrons and pos- 
itive ions is examined. 



[11] H. Ibach and H. Luth, Solid-State Physics: An Introduc- 
tion to Principles of Materials Science (Springer- Verlag, 
New York, 1995). 

[12] M. A. Ruderman, Highlights of Modern Astrophysics: 
Concepts and Controversies (John Wiley & Sons, New 
York, 1986). 

[13] V. Canuto, in Proc. 65th Enrico Fermi School: Physics 

and Astrophysics of Neutron Stars and Black Holes 

(North-Holland, Amsterdam, 1975). 
[14] E. Fradkin, Field Theories of Condensed Matter Systems 

(Addison- Wesley, Reading, 1991). 
[15] L. Pauling, General Chemistry (Dover Publications, New 

York, 1970). 

[16] G. Baym, H. A. Bethe, and C. Pethick, Nucl. Phys. 

A175, 225 (1971). 
[17] G. Baym, C. Pethick, and P. Sutherland, Astrophs. J. 

170, 299 (1971). 
[18] G. Borner, On the Properties of Matter in Neutron Stars, 

Springer Tracts in Modern Physics vol. 69 (Springer- 

Verlag, Berlin, 1973). 



* Email: gelb@alum.mit.edu. 
^ Email: kennedy@phys.ufl.edu. 

[1] J. M. Gelb, K. S. Gopinath, and D. C. Kennedy, Univ. 
of Flori da preprints UF-H EP-IFT-97-2/6; Los Alamos 



astro-ph/9702014; submitted to Phys. Rev. D 



archive 
(paper I) 

[2] S. Chakrabarty, D. Bandyopadhyay, and S. Pal, Phys. 

Rev. Lett. 78, 2898 (1997). 
[3] N. K. Glendenning, Compact Stars: Nuclear Physics, 

Particle Physics, and General Relativity (Springer- 

Verlag, New York, 1997). 
[4] S. L. Shapiro and S. A. Teukolsky, Black Holes, White 

Dwarfs, and Neutron Stars: The Physics of Compact Ob- 
jects (John Wiley & Sons, New York, 1983). 
[5] F. C. Michel, Theory of Neutron Star Magnetospheres 

(University of Chicago Press, Chicago, 1991). 
[6] L. D. Landau and E. M. Lifshitz, Quantum Mechanics: 

Non-Relativistic Theory (Pergamon Press, Oxford, 1977). 
[7] L. D. Landau and E. M. Lifshitz, The Classical Theory 

of Fields (Butterworth-Heinemann, Oxford, 1975) . 
[8] M. S. Longair, High Energy Astrophysics (Cambridge 

University Press, Cambridge, 1992). 
[9] P. Goldreich and W. H. Juhan, Astrophys. J. 157, 839 

(1969). 

[10] V. B. Berestetskii, E. M. Lifshitz, and L. P. Pitaevskii, 
Quantum Electrodynamics (Pergamon Press, Oxford, 
1982). 



8 



2000^ 

1500-^ 

1000^ 

500^ 


8 



1^ 

0.5- 

0- 

-0.5v 

-1 
8 



5000^ 
4000- 
3000 V 

2000^ 

1000- 


8 



0.1046> 

0.1046^ 

0.1045- 

0.1044- 

0.1044 
8 



2500^ 



2000- 




1 \ \ \ \ \ \ r \ \ \ \ \ \ \ \ \ \ n 



J \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ L 



CD 



00 

o 



d 



CO 

d 



CO 

d 



d 



d 



o 



u 



